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ABSTRACT 

In light of the recent inference of a high optical depth to Thomson scattering, r, from the WMAP 
data we investigate the effects of extended periods of partial ionization and ask if the value of r 
inferred by assuming a single sharp transition is an unbiased estimate. We construct and consider 
several representative ionization models and evaluate their signatures in the CMB. If r is estimated 
with a single sharp transition we show that there can be a significant bias in the derived value (and 
therefore a bias in erg as well). For WMAP noise levels the bias in r is smaller than the statistical 
uncertainty, but for Planck or a cosmic variance limited experiment the t bias could be much larger 
than the statistical uncertainties. This bias can be reduced in the ionization models we consider by 
fitting a slightly more complicated ionization history, such as a two-step ionization process. Assuming 
this two-step process we find the Planck satellite can simultaneously determine the initial redshift of 
reionization to ±2 and t to ±0.01 Uncertainty about the ionization history appears to provide a limit 
of ^ 0.005 on how well r can be estimated from CMB polarization data, much better than expected 
from WMAP but significantly worse than expected from cosmic- variance limits. 
Subject headings: 



1. INTRODUCTION 

Understanding the reionization of the intergalactic 
medium is important for cosmology for at least two rea- 
sons. How reionization occurred provides crucial data 
on the first (and possibly second) generation of sources, 
while Thomson scattering of the photons of the cos- 
mic microwave background (CMB) by free electrons sup- 
presses "primordial" anisotropics in the CMB imprinted 
at z ~ 1000. 

The first sources are not well understood (see Barkana 
& Loeb 2001 for a recent review), but there are sev- 
eral relevant observational constraints. The observation 
of quasars at z ^ 6 that show strong HI absorption 
(jBecker et al. ll2b0H) indicates that the universe had at 
least 1% of the t otal hydrogen content in neutral form 
(jFan et al. 112002^) at 2; ~ 6, with this neutral mass frac- 
tion rapidly decreasing at lower redshifts (Songaila & 
Cowie 2002). This is a strong indication that the epoch 
of reionization ended at z ~ 6. On the other hand, the 
observed anisotropics of the CMB indicate that the to- 
tal optical depth to Thomson scattering is not extremely 
high, suggesting that reionization couldn 't have started 
at redshifts much higher than about 30 l|Spergel et al. I 
|2QD3J. 

What happened between redshifts 6 and 30 is un- 
known. There has been extensive modeling and numeri- 
cal simulation, but without a good understanding of the 
sources robust conclusions are difficult to draw. Many 
studies have con cluded that reionization should hap- 
pen fairly rapidly ijCen and McDonald II2002I: iFan et aL\ 



120021) , but several recent studies have suggested (Wyithe 
& Loeb 2003; Gen 2003; Haiman & Holder 2003: here- 
after paper I) that reionization is an extended process, 
perhaps even with multiple epochs of reionization. 

In recent work (Kaplinghat et al. 2003; hereafter K03) 
it was shown that a two-step reionization process could 
have an observable signature in the large-angle CMB po- 
larization anisotropics that would provide unique infor- 
mation on the process of reionization in this difficult- 
to-probe redshift range. In this work, we extend these 
calculations to physically motivated reionization histo- 
ries provided by semi-analytic models and address the 
possibility that non-trivial ionization histories can intro- 
duce a bias in estimates of the optical depth r if a simple 
one-step model is used to fit the data. 

The results from WMAP have opened a new window 
on the dark ages. Several key cosmological parameters 
have been measured to high precision and WMAP has 
observed the si gnature of free elec trons at z ~ 10 for 
the first time (iK ogut et al. ' '2003^ . Previously, there 
had been practically no information about the ionization 
state of the IGM for 6.3 < z < 30. 

The optical depth to Thomson scattering is an impor- 
tant cosmological parameter. The temperature power 
spectrum only allows constraints on the normalization 
of the primordial gravitational potential power spectrum 
in the combination P^e"^^, so a determination of r 
allows a determination of the amplitude of potential (and 
mass) fluctuations. There has been much recent contro- 
versy over the amplitude of mass fluctuations on scales 
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of 8h^^ Mpc, CTs, as sum marized in recent parameter es- 
timates (jSnergel et al. I l2003'l. Estimates have varied by 
nearly a factor of two between different methods of deter- 
mination within the last few years, so precise and accu- 
rate CMB estimates would be invaluable. With accurate 
determinations of dg it should be possible to put stron g 
constraints on the nature of the dark energy llHu II2002D . 

In the next section we describe the models that are 
used to generate a zoo of ionization histories, while § 
3 demonstrates the effects of non-trivial ionization his- 
tories on CMB polarization anisotropics and parameter 
estimation. In § 4 we provide a crude example of how it 
is possible to reduce any biases induced by an unknown 
ionization history, and we close with a discussion. 

As a fiducial model we assume cosmological parameters 
= 0.29, f^A = 0.71, nth'^ = 0.024, h = 0.72 and an 
initial matter power spectrum P (k) (x fc in agreeme nt 
with recent results from WMAP ijSpergel et al. ll200^ . 

2. REIONIZATION MODELS 

We employ semi-analytical models of reionization, in 
order to derive a sample of physically motivated ioniza- 
tion histories. Details of the models are laid out in paper 
I, and we only briefly review the main framework here. 
The models assume that the total volume fraction of ion- 
ized regions is being driven by ionizing sources located in 
dark matter halos whose abundance is described by the 
N-body simulations of Jenkins et al. (2001) We distin- 
guish dark matter halos in three different ranges of virial 
temperatures, as follows: 

100 K <, Tvir < 10-^ K (Type II) 

10-* K <, T™. ^2x105K (Type la) 
Tvir ;> 2 X 105 K (Type lb) 

We will hereafter refer to these three different types of 
halos as Type II, Type la, and Type lb halos. Each type 
of halo plays a different role in the reionization history. 
In short, Type II halos can host the first ionizing sources, 
but only in the neutral regions of the IGM, and only if 
H2 molecules are present in sufficient quantity to allow 
efficient cooling; Type la halos can only form new ioniz- 
ing sources in the neutral IGM regions, but irrespective 
of the H2 abundance, and Type lb halos can form ioniz- 
ing sources regardless of the H2 abundance, and whether 
they are in the ionized or neutral phase of the IGM. 

The contribution of each halo to reionization is quan- 
tified by explicitly computing the expansion of the ion- 
ized Stromgren region, dictated by the source luminos- 
ity and the background IGM density and the clump- 
ing factor Chii. We allow the three different sources 
above to have different efficiencies e of injecting ioniz- 
ing radiation into the IGM. Here e = iV-y/esc/*, where 
= M*/(ribMhaio/i^m) is the fraction of baryons in 
the halo that turns into stars (^ 10% in normal stars 

0.01 in Type II halos); is the mean number of ion- 
izing photons produced by an atom cycled through stars, 
averaged over the initial mass function (IMF) of the stars 
4000 for a normal Salpeter IMF, and up to a factor 
of 20 higher for a population of massive, metal-free stars 
(Bromm, Kudritzki & Loeb 2001; Schaerer 2002); and 
/esc is the fraction of these ionizing photons that escapes 
into the IGM (-- 10% for Types Ia,b halos, and 1 for 
Type II halos). In our models, we also allow the pos- 
sibility that radiative feedback effects photo-dissociate 



II2 molecules below some critical redshift Zuv, and we 
self-consistently exclude the Type II and Type la ha- 
los from forming any ionizing sources inside regions that 
had already been reionized. We adopt a fixed Chii = 10 
in all models; variations in the clumping factor can be 
absorbed into changes in the efficiencies with redshift. 

In summary, our model has four parameters: the over- 
all efficiencies, en, eia, ejb, and the redshift Zuv at which 
II2 dissociative feedback sets in. We here use 5 differ- 
ent models that are broadly representative. Three of 
the models were chosen to have optical d epth equal to 
that measured by the WMAP experiment ijKogut et al. I 
I200.1D and two others were chosen to investigate the ef- 
fects of larger or smaller optical depths. Details of the 
models can be found in paper I. 

Model 1 assumes that for massive halos the stellar IMF 
is not metal-free (pop II) while minihalos (cooled by H2) 
form metal-free stars that produce 20 times more ion- 
izing photons. In terms of model parameters, this cor- 
responds to (ehii, eiajb) = (200,80). It is assumed that 
H2 starts to be destroyed at Zuv = 17. Model 2 assumes 
that minihalos do not contribute to reionization (efhcient 
destruction of H2) and that the efhciency in larger halos 
is increased to eia ^ 4800. In Model 3 it is assumed (as 
in Wyithe & Loeb 2003; Gen 2003) that there is a sharp 
transition from metal-free to normal stars at z = 14. 
Model 4 assumes that minihalos are more effective at 
forming stars than in Model 1, with en = 1400 and has 
molecules being destroyed at z^v ^ 25. Finally, Model 5 
assumes that feedback from star formation becomes effi- 
cient at destroying molecules at ^21 but the same 
efficiencies as Model 1. 

The ionization histories in our models are shown in 
Figure 1. As discussed in depth in paper I, the physics of 
reionization is rich in features that can naturally lead to 
distinctive ionization histories. These features can arise 
because of (1) the different types of coolants in halos with 
virial temperatures above and below ^ lO'^K, (2) the 
different response of different halos to radiative feedback 
on the II2 chemistry, and to photoionization feedback on 
gas infall, and (3) the different properties of metal-free 
and normal stellar populations. 

In all models, we assume singly-ionized Helium traces 
ionized hydrogen (i.e., 1.08 free electrons per hydrogen 
atom for a completely ionized universe). In what follows, 
discussion of ignores Helium (e.g., complete ionization 
is referred to as ~ 1), but the factor of 1.08 is included 
in all calculations. 

3. LARGE ANGLE CMB POLARIZATION ANISOTROPIES 

We modified CMB Fast ^ ijSeliak and Zaldarriagal 
11996(1 to use the ionization histories from the previous 
section to generate temperature (TT), polarization (EE), 
and cross anisotropy (TE) power spectra, CTt^CEi and 
Cap, re spectively. A similar modification was done by 
iBruscoli. Ferrara. and Scannap ieco,. (,2QQ2) but for an 
ionization history extracted from a num erical simula- 
tion, and bv lNaselskv and Chiang I l|2003j) using different 
ionization histories. The ionization histories are shown 
in Figure 1 and the corresponding power spectra are 
shown in Figures 2 and 3. 

Note in Figure 2 that for larger optical depths, there 

^ available at |http: / /www.cmbfast.org| 
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Fig. 1. — Ionization histories for the five models (explained in 
the text) considered. 



are secondar y bumps in t he Cp i- These are thoroughly 
explained in I Zaldarriagal l)1997ji . For the Cci (shown in 
Figure 3) an important difference is that we are correlat- 
ing quantities which, for each value of fc, have different 
angular frequencies on the sky. The polarization has an- 
gular frequency I = fc(ryo — rjri) (where rjo and Tjj'i are 
the conformal times today and at the onset of reioniza- 
tion) since it is projecting from the epoch of reionization 
where it was created. The temperature has a correspond- 
ingly higher angular frequency since it is projecting from 
the (further) last-scattering surface. The matched angu- 
lar frequencies of E correlated with E lead to secondary 
peaks in Cei, whereas the mismatched angular frequen- 
cies of T and E wash out the fluctuation power and do 
not lead to secondary peaks in Cci ■ 

As found in K03, for highly sensitive experiments ap- 
proaching the cosmic- variance limit, almost all the sen- 
sitivity to T comes from Cei- This is because Cei oc 
whereas Cci oc r and because the fractional uncertainty 
in Cei is smaller than the fractional uncertainty in Cci 
in the cosmic variance limit. These fractional uncertain- 
ties woul d be eq ual in the limit of perfect correlation 
{Cci = vCeiCti)- 

We nor malize P^{k = 0.05Mpc"^)e-2^ by requiring 
that the temperature fluctuation at £ = 150 is 150 fiK. 
This choice is arbitrary, but largely irrelevant for our pur- 
poses. Varying the ionization history with this product 
fixed produces no change in the angular power spectra 
at £ > 50 (except due to non-linear effects at ^ > 2000). 
Note that since Cti has been measured well in this range, 
a higher optical depth requires a larger normalization of 
Pcj) in order to agree with the data. As outlined in K03, 
variations in the fiducial model, such as a slight tilt or 
a slightly different normalization, will not have a large 
effect on our conclusions. 

To explore questions of bias in r we see how well purely 
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Fig. 2. — Polarization power spectra for the five models (see 
Figure 1) considered. Model 1 is shown as a bold dashed line, 
model 2 is shown as a dashed line and model 3 is the solid line 
with cosmic variance error bars shown. Note that models 1-3 all 
have the same optical depth. For models 4 and 5 the best fit single 
transition model polarization spectra are shown as dashed lines 
of the same color. All models are normalized to give the same 
temperature power spectrum at £ > 50. 
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Fig. 3. — Same five models as in Fig. 2 but now for the cross- 
correlation and with WMAP measurements plotted. As before, top 
line is model 4, bottom line is model 5, and models 1-3 are nearly 
indistinguishable. 



phenomenological models with one or two sharp transi- 
tions can be used to fit our physical models (l)-(5). For 
measurement uncertainties we assume Gaussian, white 
detector noise and ignore beam effects. We calculate the 
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likelihood of the phenomenological models, given one of 
the physical models as the "data", denoted now with a 
d superscript (see K03): 



X 



-21n£ 



{2e + i)fsky H[CEi + w-^]G 



p2 



Ci 



Et 



El ■ 



Tl ■ 



In the above, fsky is the fraction of sky coverage (which 
we take to be unity), and w is the weight per unit solid 
angle for polarization measurements. We have assumed 
that for our purposes we can ignore detector noise in 
the measurement of the temperature power spectrum. 
Note that we have included the effects of all the data, 
i.e., EE, TT, and TE correlations are all implicitly taken 
into account. We have not included possible effects of 
foregrounds. For WMAP, we assume w = 2 x lO^"*, 
roughly the expected two-year two channel sensitivity 
(allowing other frequencies to be used for foreground re- 
moval) while for Planck we assume w = 1.67 x 10^^, 
corresponding roughly to one year and two frequencies. 
We adopt these values as representative of the expected 
performance of the instruments, but should be viewed 
as order of magnitude estimates. For an ideal cosmic 
variance limited experiment — 0. We require that 
Cti > O.Olw"^ in the fiducial model to be included in 
the likelihood calculation to suppress contributions from 
points with low signal-to-noise. 

For simplicity, we only include terms up to £ ~ 50. 
There is practically no information in higher multipoles, 
although there is likely to be some signature a.t £ > 2000 
from non-linear effects, which could be important. We 
minimize this function by adjusting the transition red- 
shift, Ztr, of a model with sudden reionization and calcu- 
late the difference in of this best-fit sudden model 
relative to the true model. The true model is thus 
exp(Ax^/2) more likely than the most likely phenomeno- 
logical model. For Gaussian statistics, our estimator 
is equal to the usual statistic, so a rough estimate of 
the number of "sigmas" is y^Ax^. Large values of this 
misfit statistic indicate that the true model is a much 
better fit than the model being considered while small 
values indicate a model that is virtually indistinguishable 
from the input model. The best fitting sudden models 
arc indicated in Table 1, along with the difference in x^ 
and the optical depth of both the best fit and the input 
model. The first two columns indicate model number 
(see Figure 1) and true optical depth. Columns 3-5 show 
results of fitting a single sharp reionization assuming cos- 
mic variance error bars and indicate the best fit single 
reionization redshift, best fit optical depth, and differ- 
ence in x^ relative to true model. Columns 5, 6 and 7 
show best fit single transition redshift, optical depth and 
misfit statistic assuming WMAP noise levels, while the 
last three columns show the same parameters assuming 
Planck noise levels. 

For some of the models the misfit is very large, in 
one case a shift in the misfit statistic of more than 100 
(roughly "10 tr") for a cosmic variance limited exper- 
iment. This confirms that there is significantly more 
information in the large angle polarization signal than 
simply the optical depth, as shown in K03. For the 
most basic reionization signal the misfit is just above 



the "3 tj" level for cosmic variance limits, indicating 
that if reionization happens fairly quickly the exact na- 
ture of the transition is unimportant and the domi- 
nant effect on the CMB signal will be only that of the 
) op tical depth. This corresponds to t he ca se studied 
bv lBruscoli. Ferrara. and Scannapieco I ll2002l) . although 
even for this case our results are slightly less pessimistic. 
(1) Part of this is due to the higher optical depth of our 
fiducial model, providing more signal. 

From the point of view of parameter estimation, it is 
striking that the best fit sharp transition can provide a 
biased estimate of the optical depth, especially compared 
to the statistical uncertainties. For the assumed sensi- 
tivity of WMAP the statistical uncertainty St should be 
^ 0.02, for Planck St 0.005 and for cosmic variance 
errors St ^ 0.002 — 0.003. For the more exotic ionization 
histories the optical depth can be seen to be significantly 
biased for cosmic variance level measurements, with the 
direction and the magnitude of the bias sensitive to the 
details of the ionization history. Using the incorrect ion- 
ization history for model fitting introduces a systematic 
error in the value of r with a direction and magnitude 
that depends on the details of Xe{z). As seen in mod- 
els 3 and 4, offsets between the true and derived optical 
depths could easily be >0.01 for Planck. Determination 
of the optical depth, and thus the matter power spectrum 
amplitude, will be limited by a lack of understanding of 
the nature of the reionization process. For WMAP, it 
appears that the derived optical depth assuming a sin- 
gle sharp transition will not be highly biased for any of 
the reionization models that we consider, given that the 
estimated uncertainty in r is 0.03 (K03). 

4. TOWARD UNBIASED OPTICAL DEPTH ESTIMATES 

There is information in the shape of the large angle po- 
larization power spectra, so it is informative to see what 
can be gleaned. In Figure 4 we show the results of an 
analysis similar to that of K03. We used a simple two- 
step model, where it is assumed that a transition from 
partial ionization to full ionization occurred at z — 6.3 
and that a transition from nearly zero to an intermedi- 
ate (constant) ionization fraction, Xg, occurred at z — z^. 
We examine two fiducial models, one with z^ — 25 and 
the other with z^ = 16, with Xe values chosen so that 
both have r = 0.148. We then investigate the accuracy 
with which r and Zg can be recovered assuming noise 
levels typical for WMAP or expected for Planck by once 
again taking the fiducial model as the data and evalu- 
ating the likelihood (Eq. 1) as a function of the two 
parameters of the model. Figure 4 displays contours of 
constant likelihood for the three experimental cases. At 
noise levels appropriate for WMAP it will be very diffi- 
cult to differentiate between different models that yield 
the same optical depth (as pointed out in K03 and ver- 
ified by Kogut et al. 2003), but Planck would be able 
to determine the onset of partial reionization quite well 
and a cosmic variance limited experiment would be able 
to determine this onset very precisely. There will there- 
fore likely be suggestions in the data itself pointing to 
better models for measuring t. 

In the previous section, we saw that the largest biases 
arise in the cases where the misfit statistic is also large. 
Therefore when the measured r is highly biased the qual- 
ity of fit will probably also be bad, indicating a possi- 
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Table 1 Best fit single-step reionizations (roughly, number of a is -^Ax^) 
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Fig. 4. — Constant likelihood contours at 10% (thick curves) and 
1% (thin curves) of the maximum likelihood - which occurs at the 
fiducial models labeled by A and B denoted by asterisks. A is a 
model completely reionized below z = 18 and B has Nearly = 32 
and Xe = 0.4; both models have t = 0.148. All the panels have 
the same fiducial models. The solid curves correspond to model 
A, the dashed curves to model B. The upper panel is for a cosmic 
variance limited experiment (about 25 times the raw sensitivity of 
Planck). The middle panel assumes 1 Planck channel sensitivity 
while the lower panel is for two channel WMAP sensitivity. 



bly contaminated result. We now investigate whether a 
slightly more complicated fitting form for the ionization 
history can lead to better estimates of the optical depth. 
As a simple example of a path to a possibly less biased 
estimate of the optical depth, we fit models with a two- 
step reionization process, where the ionization history is 



characterized by a redshift of first ionization z^ariy, when 
the ionized fraction went quickly from effectively zero to 
Xe, and a second redshift, ziate when the ionized frac- 
tion went quickly to unity. To ensure stability in the 
numerical implementation we have jumps in the ioniza- 
tion fraction take place over a range in redshift of Az = 1 
centered on the nominal redshift of the transition and in- 
terpolated in log{z). For each ionization model from the 
previous section, we vary the three ionization parame- 
ters to find the values that minimize the misfit statistic, 
assuming cosmic variance errors only. 



Table 2 Best fit two-step reionizations, assuming 
COSMIC variance limited errors 
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Table 2 shows the derived optical depths from fitting 
a two-step model. Column 2 indicates the true opti- 
cal depth, the third column indicates redshift at which 
Xe = 1 and the fourth column shows redshift at which 
Xe changes from effectively zero at higher redshift to the 
value shown in column 5. The fifth column shows opti- 
cal depth of this best fit model, while the sixth column 
shows the difference in —2ln£ of this best fit relative to 
the true model. The final column shows the uncertainty 
in the determination of r when the extra parameters are 
allowed in the fit. This uncertainty was determined from 
the smallest and largest values of r found in models with 
Ax' < 1. 

As can be seen in Table 2, in this case the misfit statis- 
tics are much lower and the optical depth estimates are 
much less biased. In most cases the estimated values of r 
are now biased at levels near or below the cosmic variance 
statistical errors, indicating that iniperfec;t knowledge of 
the ionization history is not a fundamental limit to a 
good measure of r. In model 4 our two-step model may 
not be adequate, in that the bias in t is comparable to 
the statistical uncertainty. If the measured optical depth 
is very high (> 0.2) some care will be required to obtain 
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a precise and accurate estiraate of the amplitude of the 
matter power spectrum. Note that the reduction in bias 
has come with the cost of statistical errors in r increasing 
by a factor of 2. 

In general, ionization histories with widely separated 
(in redshift) episodes of ionization seem most prone to bi- 
ased estimates of r. We have not provided an exhaustive 
exploration of the possible parameter space of ionization 
histories, so it is still possible that reionization histories 
exist with larger biases and/or there are cases where a 
two-step reionization history does little to improve esti- 
mates of r. 

It is likely that there is a more physical parameteriza- 
tion of the ionization history that can both minimize bi- 
ases in parameter estimates and provide insight into the 
first generation of sources. If the optical depth is mea- 
sured to be higher than 0.1, as hinted by recent WMAP 
results, it will be important to find a good parameteri- 
zation. 

5. DISCUSSION 

We have shown that large angle polarization measure- 
ments could be very useful for shedding light on the end 
of the dark ages, a topic addressed in further detail in pa- 
per I. Conversely, it appears that at least a rudimentary 
understanding of the dark ages, beyond a simple optical 
depth to some characteristic redshift, will be required to 
be able to measure the amplitude of primordial fluctua- 
tions to very high accuracy. Particularly if CMB mea- 
surements suggest that the optical depth is high (> 0.2), 
there is a real danger of mis-estimating the true optical 
depth by a significant amount (> 0.01). In the near term, 
it appears that estimates of the optical depth based on 
the WMAP satellite will not be heavily biased at the level 
of the expected precision of 0.03 (K03). Therefore, the 
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derived values of erg in recent work are robust ( within 
the statistical uncertainties) to the choice of ionization 
history that is used to do the fit. This is unlikely to be 
the case for the next generation of instruments, given the 
apparently high optical depth measured by WMAP. 

Allowing even a moderately more complex ionization 
history allows much of the bias to be removed, at the 
expense of larger uncertainties. Uncertainty in the ion- 
ization history appears to provide a floor of ~ 0.005 on 
how well we can measure the optical depth to Thomson 
scattering from CMB polarization observations. 

Foreground contamination of large angle measure- 
ments is currently unknown for polarization experiments. 
With its extensive frequency coverage, Planck will be an 
exquisite instrument for characterizing and assessing the 
importance of astronomical sources of polarization. We 
believe we have been fairly conservative in our estimates 
of how many frequencies will be available for CMB mea- 
surement. 

CMB polarization measurements provide a unique 
complement to absorption studies. While hydrogen ab- 
sorption studies are sensitive to the fraction of neutral 
hydrogen, CMB polarization is sensitive to the fraction 
of ionized hydrogen. The signature of partial ionization 
at z ~ 15 will be very difficult to detect using absorption 
studies, so it will be extremely useful to have a comple- 
mentary tool. 
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